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This paper studied the reactive power control of the voltage source 
converters based high-voltage direct current system (VSC-HVDC). PI 
(proportional & integration) controller was used in this work to enhance the 
dynamic response of the system. Gravitational search algorithm (GSA) and 
sine cosine algorithm (SCA) are used to get optimal parameters of the PI 
controller. GSA algorithm is based on the gravity law for Newton while SCA 
depends on mathematical model based on cosine and sine functions. These 
algorithms have an efficient global Search capability. The VSC-HVDC is 
exposed to different disturbances for checking the controller robustness. First 
disturbance was applying three phase faults on the system. While the 
second one was applying a step change in AC voltage. Finally, 
Applying step change in regulators reference values. Simulation results 
proved the controller superiority also verified the enhancement of the system 
dynamic response. 

Copyright © 2019 Institute of Advanced Engineering and Science. 

All rights reserved. 


Corresponding Author: 

Mohamed S.Ghayad, 

Department Physics and Engineering Mathematics, 
Ain Shams University, 

Abdo Basha square, Abbassia, Cairoll517, Egypt. 
Email: M.elsayed@eng.asu.edu.eg 


1. INTRODUCTION 

The electrical power is generated, transmitted and distributed as an AC power. This power is 
transmitted from the generators to the customer via transmission lines. There are two kinds of transmission 
lines used for high voltages and long distance: HVDC and High Voltage Alternating Current (HVAC). 
HVDC transmission has many merits such as, it can transmit bulk power over long distances with minimum 
cost because it depend on two conductors only, it has minimum losses because there is no skin effect and it 
can connect two systems with different frequencies without the need of synchronization. Nowadays there are 
many projects use HVDC transmission especially in connecting different countries to exchange electrical 
power. There are two types of HVDC systems. 

a. Classical HVDC System: in this system the thyristor technology used to convert AC to DC (rectifier 
process) and from DC to AC (inverter process). The thyristors are connected in series construction to 
build up thyristor valve. The operating frequency for this valve is the network frequency (50 or 60) hertz 
(Hz) [1]. 

b. VSC based HVDC System: in this system the semiconductor used to the conversion process is the GTO 
(Gate Turn-Off Thyristor) or the IGBT (Insulated Gate Bipolar Transistor). The converters operate using 
Pulse Width Modulation (PWM) control technique [2]. 

In this paper VSC-based HVDC system is used due to its multiple advantages such as: 
a. The converters are able to independently control the reactive power and active power. 
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b. Ability of connecting passive networks which has no local power generation or weak AC networks. 

c. It has faster dynamic response due to increased switching frequency (PWM). 

d. The dc link voltage polarity remains constant during reversing the power flow process. 

The converters of VSC-HVDC system can be controlled through a closed loop consists of control 
units and PI controller which has strong robustness and simple structure. Therefor it is necessary to find 
suitable values of PI parameters for the steady operation of the system. Controlling of active and reactive 
power is one of the important requirements of the HVDC system. In classical HVDC system the reactive 
power is always a demand from the converter. Therefor it is supplied by sources of reactive power which 
installed at converter stations or by grid [3]. On the other hand the converter of VSC- HVDC has ability to 
supply or absorb reactive power into or from the grid so this system can control the reactive power of the AC 
systems at both ends independently while the flow of active power is balanced between these ends [4]. 

The aim of this paper is finding optimal parameters of PI using optimization methods to control the 
reactive power in the system during any disturbance in order to get better and faster dynamic response. The 
disturbances used in this paper are three phase fault, AC voltage sag and step change in regulators reference 
values. Conventional methods such as Smith and pole assignment methods are used to set PI parameters [5]. 
Due to the difficulty of getting the transfer function of VSC HVDC system, these methods will not useful. 
Jing Yong use estimation method to set gains of PI controller using PSCAD/EMTDC software, but the 
estimation method has no scientific foundation [6]. 

A simplex algorithm is presented for optimizing the parameters of PI controllers [7]. Modeling of 
the VSC- HVDC system for optimal power flow (OPF) solution using Newton-Raphson method [8]. 
Recently there are new methods for optimization called meta heuristic methods which has ability to get the 
optimal value of any parameter near global value fast. Genetic algorithm is used to get the parameters of 
controller [9]. Artificial bee colony and mine blast algorithms are applied to get optimal gain of VSC- HVDC 
system with sliding mode control [10]. In this paper the meta heuristic methods used are (GSA & SC A) for 
getting optimal values of the gains. 

The rest of paper is organized as follows: In Section 2, Modeling and description of the system. 
Control System is introduced in section 3. Optimization methods are presented in Section 4. Simulation 
results are provided in Section 5 and Section 6 includes conclusion. 


2. MODELING AND DESCRIPTION OF SYSTEM 

The VSC based HVDC transmission system between two AC systems consists of two converter 
stations (one is as a rectifier and the other is inverter) connected by a DC cable. Figure 1 shows the single 
line diagram of VSC-HVDC. In this paper VSC-HVDC system model which is established in MATLAB 
Simulink is used .This system is designed to transmit 200 MW at 100 KV DC. Figure 2 shows the Simulink 
model of the system. 



Figure 1. SLG for VSC-HVDC 
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Figure 2. Simulink model of VSC-HVDC 
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This model consists of: 

a. Two identical AC System (1 &2) and their ratings are 2000 MVA, 230 kV and 50 Hz. 

b. Transformers: the function of the transform is stepping down the ac system voltage from 230 kV to 100 
kV to be proper for each converter. 

c. AC Filters: Shunt high pass filters with 40 MVAR are installed after the transformers on the AC system 
side to eliminate high order harmonics. 

d. Phase Reactors: These reactors are connected in series between the transformer and the station converter 
and their values are 0.15 pu. They are used for controlling reactive and active power by regulating 
currents which passes through them and used to reduce the high frequency harmonics. 

e. Converters stations: There are two converters one acts as a rectifier and the other as an inverter and the 
semiconductor used is IGBT with antiparallel diode. There is a controller circuit in each station to reach 
controlled pulses to IGBT. 

f. DC Capacitors : These capacitors are installed on the DC side and the size of them depend on the required 
DC voltage .The function of these capacitors is reducing the voltage ripples on the DC side and used as 
energy storage to control the the power flow and they have value of 70 micro farad. 

g. DC Cables: The two converters stations are connected through DC cables .The power is transmitted 
through two pi sections of length 75 km. The construction of station 1 is shown in Fig 3 and it is identical 
for station 2. 

The construction of station 1 is shown in Figure 3 and it is identical for station 2 



3. CONTROL SYSTEM 

The controller design of converter 1 and converter 2 are identical. The two controllers are 
independent with no communication between them. These are used to control active power ( P), reactive 
power( Q) and DC voltage( U d ) in the system. Figure 4 shows an overview diagram of the control system and 
its connection with the main circuit. 



Figure 4. Overview diagram blocks of the VSC control system 
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The control circuit consists of main blocks as follow: 


3.1. The phase locked loop (PLL) block 

The function of this block is to calculate the transformation angle (0) which used in the (d-q 
transformation). This is done by measuring the frequency of the system and provides this angle. In steady 
state, sin(0) is inphase with the fundamental positive sequence of ( □ - component) and phase A of the point 
of common coupling voltage (Uabc). 

3.2. Outer active and reactive power and voltage loop 

The DC voltage and active power are controlled by the controlling the power angle (8) and the 
reactive power control is done by controlling the modulation index (which is the ratio between the 
fundamental output voltage of converter and the DC voltage). The instantaneous active and reactive power of 
the inverter can be written in terms of the d-q components of the voltage and the current as follows: 

P=\{ t^vd+^Vw,) (!) 

Q = \ (-"«<«, + 


If the reference of the dq-frame is chosen such that the (q-component) of the voltage is being very 
small and can be neglected (ULq ~ 0) then in (1) and (2) indicate that the reactive and the active power are 
proportional to the q and d component of the current respectively. Therefore, it is possible to control the DC 
voltage, the active power and the reactive power by controlling the current components i vd and i vq 
respectively. The active and reactive power and voltage loop contains the outer loop regulators that compute 
the reference value of the current vector (I re f_dq) which is the input to the inner current loop. Also it contains 
PI controller with a feed forward control to increase the speed response. 


3.3. Inner current loop 

This loop is used to track the reference current vector with a feed forward scheme to achieve fast 
control of the controller at disturbances and load changes.The equation after applying KYL for each phase is: 


Uf — U v = + J!i v 


(3) 


During unbalanced operation of equation [3], the voltage drop equations over the impedance (R + 
jcoL) holds for positive and negative sequence voltages and currents. The voltages drops are described by the 
differential in (4) as follow: 
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Where Z-{p) for positive sequence and in) for negative sequence. Equation (4) can be transformed 
to the a /?-frame as follow in (5): 
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Equation (5) can be further transformed to the rotating dq-frame as follow in (6) & (7). 



( 6 ) 

(7) 
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The positive and negative sequence voltages of the VSC side are obtained from the last equations as 

follow: 


^vdp — ^fdp — ^vdp T frtjLJi^qp — L“ i v ^p 

(8) 

^vdp — ^fdp — ^vdp E ("L}ivqp — L“ i v jp 

(9) 

^vdn — ^fdn — ^vdn — (wLJiyqjj — L~ i-ydn 

(10) 

^v-qn — ^fqn — ^v-qn E — L“ i-y^n 

(ID 

The mean voltages over the sample period k to k+1 are obtained by integrating (8), (9), (10) and 
(11) from (k) Ts to (k+1 )Ts and dividing by Ts (where Ts is the sampling time). 

^vdp — Ufdp — -fifvdp E L t*>L j \-qjt ~~ “OvdpC^ "1" ^ — H’dp 00} 

(12) 

^vqp — ^fqp — ^Vqp — (^Ljl^p — — Ovqp ^ El) — iy-qp 00 } 

(13) 

^vdn — Ufd n — — (**L) \-q Ti — “fivdnCk E lJ _ H-dn 

(14) 

Uvqn — Ufijn — ^vqn E — “Ovqn^ E l) _ ivqnOO} 

(15) 


Equations (16-19) are obtained by assuming constant network voltage and linear current during one 
sample period Ts. 

Uvdp(k + 1) = Ufdpl^O — ^3 + ~{ivqp(^ + 1) + hfqpO*-)} — ^r{VilpO s + 1) — ^vdpt^)} ^^3 

u vqp ( k + 1) = UfqpfJt) - * {ivqp (b + J.) + \- qj! (t)} - y {i^p (k + 1) + \- dp (k)} "{Vqp(k + 1) - i V qpf k j} ( 17 ) 

UvuJuOt + 1) = Ug| n Cft) —^{ininCt + 1) +i v JnO i )} — T _ {ivqiiCb+ 1) + ivqnCK)} - ^{ivdnCb + 1) _ IvdntX)} ^8) 

U vql |(b + 1) = U fqn C/£) - f {VqnCk + 1) + WqnO’Oj + y{ivdn&’ + 1) + Wfc)} - jrOvqnfr + 1) ~ Wqn(^)} (19 1 

The control is based on (16-19), where the voltages and currents at time (k+1) are thus equal to the 
reference values at time (k). 

3.4. DC voltage balance control 

This block is used to keep the voltage DC side of the three-level bridge balanced (i.e., equal pole 
voltages) in steady-state. It is possible to find small deviations between the pole voltages due to changes of 
converter current or due to nonlinearity on shortage of precision in the execution of the pulse width 
modulated bridge voltage. Figure 5 shows a high-level diagram of the simulink model for VSC controller. 


4. OPTIMIZATION METHODS 

The word of optimization stands for the process of getting optimal value for the parameters of any 
system to minimize or maximize its objective function. 

The optimization problems can be classified into many categories based on if they are discrete or 
continuous, mono or multi-objective, constrained or unconstrained and dynamic or static. An overview of 
some conventional methods and meta heuristic techniques which can be used in electrical power optimization 
problems is presented in [11]. 
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Figure 5. Simulink model of control circuit 


The conventional optimization algorithms do not give a suitable solution because the search space 
increases exponentially with the problem size. Meta heuristic optimization algorithms look to the 
optimization problems as black boxes [12]. This means that getting the mathematical models is not necessary 
because such optimization methods only vary the inputs (parameters) and notice the outputs of the system for 
minimizing and maximizing its outputs. These methods try to seek good solutions (near global) for the 
problem search space. These methods are motivated from evolutionary phenomena, behavior of creatures, 
physical laws and human-related concepts. Some of the famous algorithms in each of these subclasses are 
presented as follows: 

a) Behavior of creatures-based techniques: Particle Swarm Optimization (PSO) [13], Ant Colony 
Optimization [14] (ACO) and Artificial Bee Colony (ABC) algorithm [15]. 

b) Evolutionary techniques: Differential Evolution (DE) [16] and Genetic Algorithms (GA) [17]. 

c) Physics rules -based techniques: Colliding Bodies Optimization (CBO) [18], Gravitational Search 
Algorithm (GSA) [19] and Black Hole (BH) [20]. 

d) Human-related techniques: Mine Blast Algorithm (MBA) [21], League Championship Algorithm (LCA) 
[22] and Teaching-Learning-Based Optimization (TLBO) [23]. 

In this paper a two meta heuristic search algorithms will be used which has less computational time, 
relative high accuracy and less parameters requirements. The first method is called Gravitational Search 
Algorithm (GSA) and the second is called Sine Cosine Algorithm (SCA). The two methods are explained 
briefly in the following sections: 

4.1. Gravitational search algorithm (GSA) 

This algorithm is depending on gravity law for Newton which says “The attraction force between 
the particles is inversely proportional to the square of the distance between them and directly proportional to 
the product of their masses”. 

In this algorithm, agents are considered as particles and their performance is measured by their 
masses i.e. the better solution is corresponding to the heavy masses. All these particles attract each other by 
the attraction force and this force causes a movement of all particles towards the particles with heavier 
masses while heavy masses move more slowly than lighter ones. As the time passes, all masses will be 
attracted by the heaviest mass. This mass will present an optimum solution in the search space. The principle 
of GSA is explained in details in [19]. Figure 6 shows the steps of operation in each iteration as a flow chart. 
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Figure 6. Flow chart of GSA 


4.2. Sine cosine algorithm: 

The optimization process in this algorithm starts the by a set of random solution then it begans to 
update the position of this solution. The updating position is done according to: 

X!^ 1 = X! + x nm(r z ) x |r 3 ^ - X| | (20) 

X} +1 = X[ + r x X cos(r 2 ) X |r 3 Pi* -Xf| (21) 

Where Xj t+1 is the position of new solution, X; 1 is the position of current solution in t-iteration and rj, 
r 2 , r 3 are random numbers [24]. The switching between the two equations is based on random variable (r 4 ). 
The recent equations are the reason behind the name of this algorithm. As the iteration counter increases, the 

ranges of cosine and sine functions are updated to emphasize exploitation of the search space. The algorithm 

stops the optimization process when the iteration counter becomes higher than the maximum number of 
iterations. Figure 7 shows the flow chart of the algorithm. 
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Figure 7. Flow chart of SC A 
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5. SIMULATION RESULTS 

The objective of this section is checking the dynamic reactive power response of the HVDC system 
and trying to improve this response via tuning the PI controller parameters using optimization methods. To 
check the dynamic response the system is exposed to different disturbances such as three phase fault, AC 
voltage sag and step change in reference values. Then observing the reactive power response with respect to 
the reference in each case. This section consists of two parts, The first part explain type of parameters and the 
values of the default and optimal parameters after applying optimization methods, while the second part 
conclude the comparison between the reactive power dynamic response of default system and optimized 
system. 

Part 1: Getting optimal parameters using GSA and SC A method 

In this part the GSA and SCA are used to tune the reactive power PI controller when the system is 
exposed to severe three phase fault at system. There are six control parameters (K T _regulator, k P _control and 
k!_contol for station 1 and K T _regulator, k P _control and k!_contol for station 2) which used to control reactive 
power. K^regulator used to keep the value of measured reactive power with the desired range. k P _control 
and k!_contol are used to settle the value of the measured reactive power fast. The objective function is 
minimizing the error between the measured and reference value of reactive power. Both methods will try 
randomly to get minimum objective regardless the reactive power curve higher or lower than the reference. 
The range of these parameters is obtained by trial and error. For each method, it began to assume random 
solutions (values for the six paramters) within the given range, hence the system runs with these values. Then 
the fitness function (error between the measured reactive power value and refrence one) for each solution is 
computed. The minimum value of the fitness function indicates that these solutions can be best solutions. 
Using these solutions, the method began to generate new solutions and the previous cycle is repeated again 
until the no of iteration exceeds the maximum no of iteration (200 iteration). Range, default and optimal 
values are listed in Table 1. 


Table 1. Values of Parameters 



Parameters 

Range 

Default 

GSA 

SCA 

Station 1 

Ki_regulator 

10-30 

20 

10.8334 

10.1433 


K P _control 

1-20 

3 

2.5336 

1.0326 


Ki_control 

1-20 

3 

5.2901 

3.1119 

Station 2 

Ki_regulator 

10-30 

20 

13.6473 

21.6743 


K P _control 

1-20 

3 

6.8214 

1.6100 


Ki_control 

1-20 

3 

17.0372 

9.4688 


Part 2: Dynamic response comparison 

In this part a comparison between the default system and the optimized system in case of reactive 
power response is carried out. The system is exposed to three different disturbances as cases studies as 
follow: 


Case 1: Applying three phase faults at the bus of AC system 2 at 1.5 sec and clearing it after 0.25 
sec. The reference value of reactive power (Q r ) is -0.1 pu. The reference value of active power is ramp value 
and reaches to 1 pu after one second. The refrence value of DC voltage is 1 pu. Figure 8 shows the dynamic 
measured reactive power response (Q m ) of the default system of system 2 which illustrates that the system 
takes long time to reach the reference value of reactive power after clearing the fault. Figure 9 shows the 
optimized system response using GSA which shows that the system reaches the reference value faster than 
the case of default system. Figure 10 shows the optimized system response using SCA which also has settling 
time less than the case of default system. Figure 11 shows the behaviour of all cases in one graph to get 
accurate comparison which confirms the result mentioned before that the settling time is improved and the 
settling time of GSA and SCA is almost 2.5 s and 2.9 s respectively while for default system is greater than 
three seconds. Also, it is noted that deviation from the reference value in SCA is smaller than deviation from 
the reference value in case of GSA. Figure 12 shows the response of active power of all cases in one graph. 
The response of default and GSA case and default case are similar but the response of SCA response is better 
than them because it has small devation from the reference. Figure 13 shows the dynamic response of DC 
voltage for all cases in one graph . It is noted that the default and GSA cases are better in response than SCA 
case. 


Int J Appl Power Eng, Vol. 8, No. 2, August 2019: 145 - 158 




Int J Appl Power Eng 


ISSN: 2252-8792 


□ 153 


0.5 



0.5 

0 

Us 

O .f 
•IS 


1 - T ' '1 - 

^ 1 

_ 1 _!_ 1 _ 

r ii: 

_!_i_ 

II 

cc 

■ * 


OS 


IS 

Tijurn) 


2.5 


Figure 9. Reactive power response for GSA case 
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Figure 11. Reactive power response for all cases 



Figure 12. Active power response for all cases 
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Figure 13. DC voltage response for all cases 


Case 2: Applying a step change in AC voltage of system 1 and this is can be done by using three-Phase 
programmable voltage Source in system 1. The amplitude of AC voltage of system 1 is changed from 1 pu to 
0.8 pu at 1.5 s for 0.25 s. Figure 14 shows the dynamic reactive power response of the default system of 
system 1 which illustrates that the deviation from the reference value is high during the voltage sag which 
reachs to -0.48 pu. Figure 15 shows the optimized system response using GSA which has less deviation from 
the reference value than the default case during the voltage sag and this value reach to -0.4 pu. Figure 16 
shows the optimized system response using SC A which has minimum deviation from the reference value 
compared with default and GSA cases during the voltage sag and its value reach -0.22. Figure 17 shows the 
all previous cases in one graph which confirms the result mentioned before that SCA case has less reactive 
power deviation from the reference value. Figure 18 shows the dynamic response of active power which 
indicate that the response in case of SCA has less deviation from the reference value. Figure 19 shows the 
dynamic response of DC voltage which indicate tha the AC voltage sag has small effect on DC voltage and 
the response of all cases is similar. 
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Figure 15. Reactive power response for GSA case 
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Figure 17. Reactive power response for all cases 



Figure 18. Active power response for all cases 


VP . 1 

n 



T I 



r 

1 



* 






J_i 

i 

1 1. 1 


-RcTeiTHi-r 
-I>t fa till 
-G! SA 
'S(A 


0.5 




2.5 


Figure 19. DC voltage response for all cases 


Case 3: Applying step change in regulators reference values of active power, reactive power and DC voltage. 
After the system is reached to the steady state, the DC voltage reference is changed from 1 pu to 0.85 pu at 
1.3 s then a step change in reactive power reference is applied at station 1 from 0 pu to -0.1 pu at 1.7 s last, a 
step change in active power reference is applied at station 1 from 1 pu to 0.9 pu at 2.3 s. Figure 20 shows the 
dynamic reactive power response of the default system in station 1 which illustrates that the measured 
reactive power has high deviation from the reference value. Figure 21 shows the optimized system response 
has almost no deviation from the reference value of reactive power and it reaches the reference value fast. 
Figure 22 shows the optimized system response using SC A which also has very small deviation from the 
reference value but the settling time is greater than the case of GSA. Figure 23 shows the previous cases in 
one graph which confirms the result mentioned before that both GSA and SCA can improve the dynamic 
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response during changing of regulators reference values. Figure 24 shows the dynamic response of active 
power which is similar in all cases but Figure 25 shows the dynamic response of DC voltage which is best in 
case of SCA and GSA as it reach to the reference value faster than the default case. 



Figure 20. Reactive power response for default case 
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Figure 21. Reactive power response for GSA case 
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Figure 22. Reactive power response for SCA case 



Figure 23. Reactive power response for all cases 
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Figure 24. Active power response for all cases 
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Figure 25. DC voltage response for all cases 


6. CONCLUSION 

The main purpose of this paper is improving the dynamic response of VSC HVDC transmission 
system during existing of any disturbance in it such as three phase fault, voltage sag of any AC grid and 
change in reference values. This is done by getting optimal gains valus of PI controllers to control reactive 
power. The optimization methods used to get optimal gains are GSA and SCA which Meta heuristic 
algorithms. Results proved that PI-SCA has better response than PI-GSA regarding reactive power deviation, 
but PI-GSA is has better performance in settling time. 
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